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Abstract 

It is shown that the rates of the decays Y(l 3 Z?i) —* rj T(15) and T(1 3 I?2) — * V T(15) 
should be comparable to and likely exceed that of the recently discussed in the literature 
two-pion transition T(l-D) — * it ttT(1S). The reason for this behavior is that the 
discussed rj transitions are enhanced by the contribution of the anomaly in the flavor 
singlet axial current in QCD. 



1 Introduction 



The D-wave states in the family of the bb resonances present a new interesting testing 
ground for the study of heavy quark dynamics. The 3 Dj states with J = 1,2,3 should 
form a closely spaced triplet of resonances, of which one (most likely the 3 D2) with the 
mass of about 10.16 GeV has been recently observed^ in the CLEO experiment through 
the radiative transitions to and from the D wave state. It is clear however that similarly 
to other excited bb resonances there should also be strong-interaction transitions from the 
D states to lower resonances with emission of light mesons, i.e. of two pions, rj, and also 
a weaker isospin violating transition with emission of tt°. In particular the transitions of 
the type T(LD) — ► 7T7tT(1S') were discussed in the literature [2J El H] m some detail. The 
interest to the two-pion transitions is explained by that these are well known to be dominant 
for the more familiar cases of the hadronic transitions between 3 Si states, i.e. from ip' in 
charmonium and from T(2S), while the rate of similar decays with the rj emission being 
considerably smaller, and the related isospin violating transition with emission of a single 
7T° being greatly suppressed further. The amplitudes of these transitions between the 3 Si 
states of heavy quarkonium and the pattern of the rates were understood long ago0 El Ej 
within the general method of describing the hadronic transitions in heavy quarkonium using 
the multipole expansion [SI El i n QCD for the interaction of the heavy quarkonium with soft 
gluonic field. In the transitions between the 3 S\ heavy resonances the relevant amplitudes 
for production of the light mesons are determined [H] by the low energy theorems arising from 
the quantum anomalies in QCD: the emission of an S wave pair of pions is dominated by 
the anomaly in the trace of the energy-momentum tensor ^0], while the P wave emission 
of the rj is regulated by the anomaly in the flavor singlet axial current [TT| . The anomalous 
contribution greatly enhances both rates [5] in agreement with the available data. 

The purpose of the present paper is to point out that the relation between the rates 
of the two-pion and 77 transitions from T(l 3 Di 2) to T(15) should be quite different from 
the pattern observed in the transitions between the 3 Si states. Namely, in the ID — > IS 
transitions the pions are emitted in the D wave, and the corresponding amplitude for the 
production of the pions by the relevant gluonic operator decouples from the anomaly. On 
the other hand, the P wave emission of rj is still possible for transitions to the 3 Si state from 
the l 3 Dj resonances with J = 1 and J = 2, and, as will be shown here, is indeed contributed 
by the axial anomaly in QCD. The resulting enhancement of the amplitude of the rj emission 
turns out to be sufficient to compensate for the suppression factors inherent in this decay (the 
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flavor SU(3) violation, as well as a suppression by the inverse of the b quark mass), so that 
the rate of the r\ transitions should be comparable to that of the two-pion ones, and in fact 
is quite likely to be the largest among the hadronic transitions from the 3 Dj states 1 . A more 
definite quantitative estimate of the ratio of the rates, r(l 3 _D lj2 — > rj l 3 Si)/T(lD — > n tt IS), 
is hindered by the present poor understanding of a parameter governing the non-anomalous 
amplitude of production of the pion pair by gluonic operators, as will be discussed in Sect. 4. 

Within the QCD multipole expansion treatment of the hadronic transitions in a heavy 
quarkonium, outlined below in Sect. 2, the evaluation of the amplitude of the rj emission 
requires knowledge of the matrix element (r]\G pu D p G\ a \0), where is the gluonic field 
tensor, and D p is the covariant derivative. It will be shown in Sect. 3 that this matrix 
element is completely determined by the well known expression 



following from the anomaly in the divergence of the flavor singlet axial current in QCD. In 
eq.(HJ) fn is the r] 'decay constant', equal to the pion decay constant f n ~ 130 Me V in the 
limit of exact flavor SU(3) symmetry, and f n is likely to be larger due to effects of the SU(3) 
violation. Also throughout this paper the normalization of the gluon field tensor includes the 
QCD coupling g (so that e.g. the Lagrangian for the gluon field reads as L g = — G 2 /(4g 2 )). 

The amplitudes and probabilities of specific decays are calculated in Sect. 4. Besides a 
numerical comparison of the rates for the two-pion and r\ transitions between the 3 Dj and the 
3 Si resonances in the bb system, also discussed there are the greatly suppressed transitions 
with emission of a 7r°. 

Finally, the Section 5 contains a summary of the discussion in the present paper. 

2 Transition amplitudes in the multipole expansion 

We start with a brief reminder of the leading terms in the multipole expansion in QCD which 

are relevant to the discussed transitions [SJ El- 

1 This behavior is reminiscent of that expected|12| for transitions between 1 P\ and 3 Si states. There the 
two-pion emission, also decoupled from the anomaly, is additionally kinematically suppressed, so that the 
isospin violating single 7r° emission becomes more probable due to the contribution of the axial anomaly 
(while the r\ emission is impossible kinematically). 




(1) 
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The two-pion transition arises in the second order in the El interaction with the chro- 
moelectric gluon field E a described by the Hamiltonian 

H E1 = ~Cr-E a (0) , (2) 

where £ a = — 1\ is the difference of the color generators acting on the quark and antiquark 
(e.g. t\ = X a /2 with A a being the Gell-Mann matrices), and r is the vector for relative 
position of the quark and the antiquark. 

The transitions of the type 3 Dj — > r] 3 Si are induced by the interference of the El 
interaction in eq.(j2J with the M2 term containing the chromomagnetic field B a and described 
by the Hamiltonian 

H M 2 = -(4m Q )- 1 r Sjn (A^(0)) a , (3) 

where D is the QCD covariant derivative, rnq is the heavy quark mass, and S = {B\ + d^jl 
is the operator of the total spin of the quark-antiquark pair. It should be noted that the Ml 
term, formally of a lower order in the multipole expansion, is proportional to the spin-flip 
operator {a\ — B2) and thus does not contribute to transitions between states with the same 
total spin. 

Using the expressions (J2J) and the transition amplitudes are found in the standard 
way: 

A** = A( 3 D,J ► 7T 7T ^Sl) = (7T7T|^|0) Aij , (4) 

AW = A( 3 Dj -^ V 3 S 1 ) = m Q l ( V \E? (D 3 B k y + {D s B k )* E? \ 0) , (5) 
where Ay and A^\ are the heavy quarkonium amplitudes, defined as 

A^^lSrr^r^lD) (6) 

and 

42 = li ('SilCngrjeS^Dj) (7) 

with Q being the Green's function of the heavy quark pair in a color octet state, and also 
in these expressions a use is made of the fact that in matrix elements between color singlet 
states one can replace £ a . . . £ b by (5 ab /8) £ c . . . £ c . 

The (chromo)electric dipole interaction in eq.(J2J) does not involve the spin of the heavy 
quarks. Thus in the leading nonrelativistic limit, assumed throughout this paper, where the 
spin and coordinate degrees of freedom can be considered as independent, the amplitude A^ 
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in eq.(jl|) does not depend on the spin variables of the initial D wave state or of the final S 
wave one. For this reason the rates of transitions from each of the 3 Dj states to tttc 3 Si (and 
also of 1 D 2 - * tttt 1 5'o) are all the same |1 (modulo small differences in the energy release, 
whose effect is formally beyond the assumed approximation) and can in fact be calculated 
for spinless quarks. 

The amplitude does depend on the spin-orbital state of the quark pair and is different 
for different values of the total momentum J. For the purpose of the present discussion in 
the leading approximation of the decoupled spin variables it is convenient to first represent 
this amplitude not in the basis of states with definite J, but rather in a form with explicitly 
factorized spin and orbital components in the Cartesian coordinates. For this representation 
we denote as Q and \% the spin polarization amplitude of respectively the initial 3 D state 
and the final 3 Si state, and as ip^ the orbital polarization amplitude of the L = 2 wave in 
the initial D states. The tensor ipij is symmetric and traceless, as appropriate for an L = 2 
state, and all these amplitudes are assumed to be normalized in the standard way, so that 
the sums over the polarization states are defined as 

£ x*Xj = Yl CO = Sij, Yl ^ij^ ki = 2 ( 6ik6jl + 6il6jk ~ \ 6ij6kl ) ■ ( 8 ) 

pol pol pol 

In this notation the amplitude Ay (for spinless quarks) is proportional to fyj and can thus 
be written in terms of a scalar quantity A 2 as = ipij A 2 , while the amplitude is 
expressed in terms of the same A 2 as 

42 = ~ e k imX*P {J) Cm A 2 , (9) 

where is the projector on states with definite J, acting on the product of the spin and 
orbital polarization amplitudes ipijCm- 

The quantity A 2 depends on details of dynamics of heavy quarkonium, and at present is 
highly model-dependent. For this reason a prediction of the absolute rates of the discussed 
decays involves a considerable uncertainty. Clearly, however, A 2 cancels in the considered 
here ratio of the rates of the two-pion and rj transitions, which is thus determined by the 
ratio of the matrix elements entering the equations (jlj and (|5)L describing the production 
by the gluon operators of the corresponding light meson states. 
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3 Matrix elements for production of light mesons by 
gluonic operators 



The gluonic matrix element for the two-pion production in eq.(@} multiplies the traceless 
tensor ip^ and thus receives no contribution from the (enhanced) trace anomaly in QCD. 
Rather this matrix element is parameterized [SI 12] i n terms of the QCD coupling a s and the 
parameter pa introduced in Ref.[6 j as 'the fraction of the pion momentum carried by gluons'. 
Using this parameterization, one can write 

A n+n - = (7r + 7r|£?££|0) Aj (xlCk) A 2 = Ana^optpjAj (xlCk) A 2 , (10) 

where the final state with charged pions is assumed for definiteness, and p ± stand for the 
momenta of the pions in the heavy quarkonium rest frame. 

The matrix element of the gluonic operators in eq.(jHJ) can be found as soon as the am- 
plitudes of general Lorentz structure (r]\G pu (D p G Ao-) a |0) and {r]\(D p G pu ) a G Xa \0) are known. 
These structures can in fact be reduced to (a total derivative of) the amplitude described by 
eq.((TJ). The possibility of the reduction of the structures with derivatives to the expression in 
eq-O is determined by the general theory [Hj. Here we present the specific implementation 
of this reduction, which uses the following simple algebraic identity valid for an arbitrary 
four-vector p: 

Pp^-pvXu PX^pupa Pa^pvpX Pp^upXa Pu^ppXcr j (H) 

where the convention €0123 = 1 is assumed. The antisymmetry of the field tensor G pv then 
allows one to write the general form of the first of the discussed matrix elements in the linear 
order in the rj momentum p in terms of two scalars X and Y: 

Pa^-pvpX ) • (12) 

The third structure, allowed by the symmetry and proportional to {j> p t vp \ a — p u e pp \ a ), is 
reduced to the first two due to the identity (fTTjl . Furthermore, applying in eq. ffT2|) the 
equations of motion (the Jacobi identity): D p G\ a + D a G p \ + D\G ap = 0, one arrives at the 
relation X = 2Y. 

Likewise, writing the second of the discussed matrix elements in terms of two scalars X 
and Y as 

i {rj(p)\(D p G pi/ ) a G a Xa \0} = Xp p t puXu + Y {p p t Xapv - p v t Xapp ) , (13) 
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and applying the Jacobi identity, one finds the relation X = 2Y. 

The sum of the expressions (J 12)) and (|13|) should combine into a total derivative, i.e. 
the sum should be proportional to p p . This is possible due to the identity (fTTj) under the 
condition that Y = Y, so that all the considered scalar form factors are expressed in terms 
of one of them, e.g. in terms of X: 2 

X = X, Y = Y = -X. (14) 

Using this relation and contracting the sum of the expressions (JED and (UH) with \e^ vXu the 
form factor X is identified from the equation (0) as 

1 ~ 4vr 2 [2 

X = ~30 ^\G a G a \0) = -— ^- /X • (15) 

The relations (fT2"|) - fj!4|) fully define the gluonic matrix element in eq.© in terms of X: 

1 (77 \E°; (DjB k ) a + (DjB k ) a Ef\ 0) =-2X(3p j 6 ik -p i 8 jk ) . (16) 

One can now use this expression and the form of the amplitude from eq. to write the 
full amplitude of the rj transition as 

4^ = X P l £jlmX*lP {J) Aj Cm A 2 , (17) 

where the symmetry of the tensor amplitude ipij is taken into account. 



4 Relations between the decay rates 

It is quite clear from the equation (|17|) that due to the presence of 6ji m the amplitude of the 
r) emission vanishes for the J = 3 state, which is totally symmetric in the indices i,j, I, when 
expressed in terms of the product tpijCm- This previously mentioned behavior is naturally 
expected, given that the 77 meson is emitted in the P wave. The projection of the latter 
product on the state with J = 1 is given by 

(!) ^ ^ 

2 An alternative derivation of two of these relations, namely X — X and Y — Y, would be by arguing that 
in the particular amplitudes (|12|l and (|13|l the operators G a and (DG) a can be considered as commuting 
with each other, so that the expressions (|12|l and differ only by re-labeling the indices. 
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The normalization factor here is readily found from eq.© and the condition that the sum 
over all the polarization states should be equal to the number of the polarization states of 
3 ^i: E^,jP (1 %-Cm| 2 = 3. 

Using the explicit expression (|18|) in the formula (|17|) for the transition amplitude, it is 
a straightforward exercise to find the square of the amplitude averaged over polarizations of 
the initial 3 D\ state and summed over the polarizations of the final 3 Si state: 

10 n X 2 



\ a0-)\2 — ri 2 HI 2 (-\Q) 

IA, | - g \A 2 \ , (19) 

where the overline in the l.h.s. denotes the prescribed averaging-summation operation, and 
Pv = \P\- 

In order to find the similar quantity for the transitions from the 3 D 2 states, one generally 
would have to consider the projector P^ similarly to the previous treatment of the projector 
However, at this point it is simpler to use the fact that only the J = 1 and J = 2 states 
contribute to the 'grand total' sum of the square of the amplitude over all orbital and spin 
polarization states, and to find the J = 2 contribution by using eq.([T9"]) and subtracting the 
sum over the J = 1 states from the 'grand total': 

1 /40 10\ 2 X 2 , A l2 2 X 2 
5 V 3 3 / 'rrin rnj 




\A?? = -A E KSl'-'^N ?-T P ^|^ = 2 P ^|^. (20) 



A comparison of the latter result with that in eq. immediately leads to the prediction 
of the ratio of the decay rates: 

r(3 jDl ^ v 3 5l ) = 5 r(3jD2 ^ v s Si) (21) 
y 

It should be noted that this relation is obtained in the limit where all effects of the spin- 
dependent interaction in the heavy quarkonium are neglected. The ignored effects include 
in particular the fine-structure splitting between the masses of the 3 D 2 and 3 Di. However 
this splitting is expected to be quite small, not larger than about 10 MeV (for a summary 
of potential model predictions see e.g. Ref.|15j). Thus if the difference in the kinematical 
factors p 3 in the decay rate is used as a representative measure of the contribution of the 
unaccounted corrections, one might expect that the accuracy of the relation (j2"Tj) should be 
about 10%. 

In order to estimate the relative rate of the discussed r] and tctt transitions we need to 
compare the corresponding phase space integrals of the corresponding amplitude squared at 



7 



the energy release A = M(l 3 Dj) — M(T) ~ 700 MeV, i.e. to compare the expression 
(numerically, p v ~ 435 MeV) for the rj emission with the integral 



/WMA-,-^ ^ ^ (23) 



for the two-pion emission, where in each of the integrals £ stands for the energy of the 
corresponding light meson. 

The square of the 7r + 7r~ transition amplitude from eq. (jlO|) . averaged over the initial and 
summed over the final polarization states, is given by 



I^TT+TT-I 2 = -= (a s p G ) 2 
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p\pI + ^{pi-nf 



\A 2 \ 2 ^fn 2 (a sPG y P lpj\A 2 \\ (24) 



where are the momenta of the pions, and pi j2 = \Pi,z\- Also in the last transition in eq. (|24j) 
the averaging over the relative angle between the momenta is performed, as appropriate for 
the purpose of calculating the integral in eq. (|2l?|) . A straightforward integration in eq.(J2H|) 
of the expression (|24j) yields 

WW = 0.44 {a sPG f\ A 2 \ 2 , (25) 

where the numerical factor 0.44 accounts for the finite mass of the pions in the phase space 
integral (i.e. for massless pions this factor would be equal to one). Thus using in eq.()22|l the 
expression for the rj emission amplitude from eq.([20|). the ratio of the decay rates is found as 

T(T(1 3 D 2 ) ^ V T) W® 



T(T(1 3 D 2 ) -> tt+tt-T) W, 



X 2 pl 448 vr 4 fl mtpl /0.64\ . . 

0AA(a s p G y^ml 15 0AA(a s p G ) 2 rn 2 A 7 \a s p G J ' 1 ) 

Here the expression for X given by eq. (|15J) is used, and also in the final numerical calculation 
rather conservative values of f v and m& are assumed: f v ^f n ^i 130 MeV, mi, ps 5 GeV. 

Clearly, the main uncertainty in evaluation of the ratio of the decay rates comes from 
poor knowledge of the dimensionless parameter a s p G with both factors normalized at a scale 
p, set by the characteristic size of the quarkonium. The estimates of the relevant value of 
this parameter range from a s p G ~ 0.2 [2 J to a s p G ~ 0.59 [3], with a realistic value likely 
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being close to 0.35. In either case, the numerical result in the equation ()26|) predicts that 
the r] transition rate should be not smaller, but most plausibly larger, than the rate of the 
transition with the emission of two pions. 

The presented analysis can also be readily applied for an estimate of the rate of the isospin 
violating transition with emission of a single 7r°. In the discussed axial-anomaly-dominated 
processes the amplitude of the 7r° transition is simply related?] to the amplitude of the 7] 
emission by the ratio ^T] of the corresponding gluonic matrix elements: 

A n _ (iT \G a G a \0} p n _^ rn d - m u f n m\ p n 
A v (r]\G a G a \0) p v m d + m u f v p v 

where m u and m d are the masses of the up and down quarks, and the ratio of the decay 
rates is thus estimated as 

r(T(l 3 D J )^7r T) o (m d -m u y ft m% p\ 



r(T(l3 J D J ) -^T) \m d + mj p* ' 

For (m d — m u )/(m d + m u ) ~ 0.3 this estimate gives numerically about 4 x 10~ 3 , so that 
the rate of the transition with the emission of a single pion should be quite small. 



5 Summary 

Within the description of hadronic transitions among the levels of a heavy quarkonium, based 
on the multipole expansion in QCD, the interference of the El interaction in eq.(J2J) with the 
M2 term in eq.® gives rise to transitions between the 3 Dj and 3 Si resonances in the bb 
system with the emission of the rj meson in P wave. These processes are observable as the 
decays T(l 3 Dj) — > rjT of the D wave resonances with J = 1 and J = 2. The dependence 
on the coordinates of the heavy quarks of the transition amplitude describing these decays 
coincides with that of the amplitude for the transitions T(l 3 Dj) — > 7r7rT. Thus the highly 
model-dependent factor related to the dynamics of the heavy quarkonium cancels in the 
ratio of the amplitudes, which ratio is thus determined by the amplitudes of production 
of the corresponding light meson states by gluonic operators. The relevant amplitude for 
the production of rj (eqs. (fTK| HH )l ) is shown to be completely determined by the low energy 
theorem in eq.Q directly related to the anomaly in the flavor singlet axial current in QCD. 
The enhancement of the rj emission due to the anomaly contribution, which emission normally 
would be suppressed by the inverse of the b quark mass, and the relative smallness of the 
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flavor SU(3) breaking, makes the rate of this process comparable to, or most likely greater 
than that of the two-pion emission, as described by the equation (J26j) . This pattern of the 
relative rates is qualitatively different from that observed in the 2 3 S*i — > l 3 Si transitions, 
both in charmonium and in the bb system, where the two-pion transitions dominate. The 
difference is that in the latter transitions the two-pion emission amplitude is also enhanced 
by an anomaly, in this case the conformal anomaly in QCD, while the amplitudes of the 
D — > mi S decays receive no enhanced anomalous contribution. The presented treatment, 
based on the multipole expansion, predicts the ratio of the rj transition rates to T from 
the 3 L>i and 3 D 2 resonances (eq. (|21|) ) with plausibly an accuracy of about 10%, while the 
numerical value of the ratio of these rates to that of the two-pion emission, given by eq. (|26j) . 
contains a considerably larger uncertainty associated with the present poor knowledge of the 
parameter a s pc, describing the emission of two pions (eq.(|IDj)). 
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